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ABSTRACT 

We propose a dynamical mechanism for capturing stars around a massive black hole (MBH), which is based 
on the accumulation of a very dense cluster of compact stellar remnants near the MBH. This study is motivated 
by the presence of ~ 10 young massive stars (Mi, ~ 3-15 Mq, spectral types ^B9V-08V) less than 0.04 pc 
from the MBH in the Galactic center (GC). Their existence in the extreme environment so close to an MBH 
is a challenge for theories of star formation and stellar dynamics. We show that young stars, which formed 
far from the MBH and were then scattered into eccentric orbits, repeatedly cross a cluster of stellar black 
holes (SBHs), where they may undergo rare direct three-body exchanges with an MBH-SBH "binary". The 
interaction between two objects of comparable mass ejects the SBH and captures the star on a tight orbit around 
the MBH. Such captures can naturally explain some trends observed in the orbits of the young stars. We derive 
the capture cross-section, validate it by Monte Carlo simulations, and calculate the number of captured stars in 
the GC using the currently uncertain estimates of the numbers of SBHs in the inner 0.04 pc and of young stars 
in the inner few parsecs of the GC. We find that under favorable conditions three-body exchange can account 
for ~ 25% of the observed stars, mostly at the fainter end of the observed range. We discuss additional effects 
that possibly increase the capture efficiency. Future detections of the dark mass around the MBH and deeper 
surveys of the central parsecs will establish whether or not there are enough SBHs and young stars there for 
exchange captures to singly account for the central young stars. We estimate that there are also ^ 35 lower 
mass stars (A/^, ~ 1-3 Mq, ~G2V-A0V) in the inner 0.04 pc similarly captured by exchanges with neutron 
stars (NSs). Ongoing replacement of compact remnants by main-sequence stars (SBHs by NS progenitors, NSs 
by white dwarf progenitors) may regulate the accumulation of compact remnants near the MBH. 
Subject headings: black hole physics — galaxies: nuclei — stars: kinematics 



1. EVITRODUCTION 

Deep near-infrared photometric (Krabbe et al. '1995^, Gen- 
zel et al. 2003), spectroscopic (Genzel et al. 1997; Eckart, 
Ott & Genzel 1999; Figer et al. 2000; Geza ri et a l. 120021: 
Ghez et al. 2003) and astrometric (Ghez et al. I2003t Schodel 
et al. ^003) observations of the dense stellar cusp around the 
massive black hole (MBH) in the Galactic center (GC) reveal 
a centrally concentrated distribution of young massive stars, 
iV^ ^ 10 stars within the central ^ 0.04 pc, ^ 40 within 
~ 0. 1 pc, spanning a mass range of ~ 3-15 Mq (spectral types 
-B9V-08 V) with a median mass of 10 Mq, radius of 
~ 4.5 Rq and main sequence life time of ~ 3 x 10^ yr' . 
Orbital solutions obtained for eight of the stars (Ghez et al. 
120031 Schodel et al. 2003 ) tentatively suggest some trends 
in their orbital properties: a lower bound on the apoapse of 
~ 0.01 pc (Ghez et al. 2003) and higher than random orbital 
eccentricities (Schodel et al. 2003). 

None of the solutions proposed so far for the puzzl e of 
the young stars (Genzel et al. 2003; Ghez et al. '2003"^ are 
satisfactory. These fall into three categories: exotic modes 
of star formation near the MBH; rejuvenation of old stars 
in the local population; or dynamic migration from farther 
out, where stars can form. Even if shock cooling of molec- 
ular gas by cloud-cloud collisions could trigger star forma- 

' These are rough inferences. The separation of the young stars from the 
old population is uncertain since at present only the brightest star has been 



spectroscopically identified (Ghez et al. 1200! 



tion near the MBH (Morris 119931) . molecular clouds would 
either form stars or be tidally disrupted well outside of the 
inner 0. 1 pc (Vollmer & Duschl 2001 ). Growth and rejuvena- 
tion by mergers (Genzel et al. 2003) are not expected to be 
efficient in high velocity collisions near an MBH. Tidal heat- 
ing by the MBH requires that the stars approach the MBH 
much closer than they are observed to do (Alexander & Mor- 
ris 2003). The young stars are too short-lived to have formed 
far from the MBH and then migrated inward by mass segrega- 
tion or dynamical friction. The migration can be accelerated 
if the stars are associated with a massive "anchor": an ex- 
tremely dense young cluster (Portegies Zwart, McMillan & 
Gerhard Kim & Morris 2003), a very massive binary 

companion (Gould & Quillen 2003), or a 10^-10"' Mq black 
hole (Hansen & Milosavljevic 2003). However, these scenar- 
ios must assume the existence of very rare, or even hypotheti- 
cal objects, or else they cannot bring the stars close enough to 
the MBH. Such processes may possibly explain the separate 
population of very massive and luminous "He stars" 0.1-0.5 
pc from the MBH (Krabbe et al. 1995), which we do not at- 
tempt to model here. 

Our model is based on the fact that lO'^-lO^ stellar black 
holes (SBHs) of mass ^ 7-10 Mq are estimated to exist 
within 1 pc of the MBH in the GC, where they have been 
accumulating by dynamical friction over the lifetime of the 
Galaxy (tH ^ 10 G yr ) from a "collection basin" ^ 10 pc 
wide rMoms ll993t Miralda-Escude & Gould 200 0). Nu mer- 
ical simulations of the evolution of the GC ("Freitag l2003l) con- 
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firm that the SBHs sink to the center on a short timescale of a 
few gigayears, settle into a centrally concentrated distribution 
where the enclosed number scales as iV, (< r) oc r^/** (Bahcall 
& Wolf 1977), and dominate the stellar mass there. 

2. STARS AND REMNANTS IN THE GALACTIC CENTER 

On the ~ 1 pc scale, conditions ar e mor e favorable for star 
formation. Observations (Figer et al. '1999') and theoretical ar- 
guments (Morris 1993) indicate that star formation in the GC 
is ongoing and is significantly biased toward massive stars. 
Here we represent the present-day mass function (PMF) of 
the central few parsecs of the GC by a simple (nonunique) 
model that has such properties. We assume that the PMF is the 
product of continuous star formation at a constant rate with a 
Salpeter initial mass function (IMF) in the range Mi, — 1.5- 
120 Mq, so that d7V*/dA/^(i) cx M~'^-'^^ \m-n[t,U(Mi:)\. 
Stars with ikf^ < SMq are assumed to evolve into 0.6 Mq 
white dwarfs; those with M^, = 8-30 M© into 1.4 M© neu- 
tron stars (NSs); and those with M^, > 30 M© into 1 Mq 
SBHs (consistent with the mass distribution found in black 
hole binaries; McClintock & Remillard 2004). The gas lost in 
the course of stellar evolution is assumed to be ef ficient ly ex- 
pelled from the system. Using the Schaller et al. (il992h solar 
metallicity stellar evolution tracks, we find that at t = t^f ~ 10 
Gyr, the mass fraction in stars out of the total mass is 0.22, 
the mean stellar mass is 2 Mq and the number fraction of 
young stars in the range 3-15 Mq is /j, = 0.06, with a mean 
mass of Mi, — 4 Mq (corresponding to R^, = 2.34 Rq and 
ti, = 2 X 10^ yr). Dynamical measurements of the cen- 
tral gravitational potential indicate that the MBH mass is 
m ^ 3 X 10*5 Mq (Schodel et al. 2003; Ghez et al. 2003) 
and that in the range ~ 0.5-10 pc the stellar mass distribu- 
tion is well represented by a spherically symmetric power- 
law distribution, cx r^", with the most recent calibration 
giving pi:{l pc) = 1.2 x 10^ Mq pc"^ and a = 1.7 (R. Gen- 
zel 2003, private communication). The mass distribution and 
mass function translate to 1.5x 10^ SBHs born inside the col- 
lection basin (< 5 pc) and to 2.3 x 10** young stars inside 2.5 
pc, where the enclosed dynamical mass is 3.4 x 10^ Mq, in 
general agreement with the observationally based PMF in the 
inner 15 pc of Mezger et al. ( 1999), scaled to 2.5 pc. 

Young stars formed on the > 1 pc scale will not be on ec- 
centric orbits initially, because their progenitor clouds cannot 
survive the tidal field of the MBH. Over time, the stars will 
be scattered by gravitational perturbations. The most efficient 
of these, which can operate on a timescale of less than t*, are 
due to massive star-forming clusters (Zhao, Haehnelt & Rees 
One such cluster is statistically expected to exist, un- 
detec ted, w ithin a few parsecs of the MBH (Portegies Zwart 
et al. '2002"). It is also possible, albeit speculative, that cloud- 
cloud collisions on the few parsecs scale could lead to rapid 
formation of massive stars on radial orbits. 

Genzel et al. 12003 ) estimate that the enclosed stel- 
lar mass inside r = 0.4 pc scales as Af*(< r) = 2.5 x 
10^(r/l pc)^'^^ Mq, assuming that the mass follows the star 
counts from 1 pc to the center However, the total mass dis- 
tribution (stars and compact remnants) inside 0.1 pc, where 
the MBH dominates the potential but where the remnants 
domin ate the extended mass, is not known (Mouawad et al. 
120041) . The central density of the SBH cluster depends on 
various uncertain quantities: the SBH mass function, the stel- 
lar IMF and formation rate, the remnant progenitor masses, 
and the dynamical age of the GC. A rough upper limit on 
N,{< r,t) can be obtained by requiring that the SBHs sur- 



vive being drained into the MBH over the lifetime of the sys- 
tem, dN,/dt < N,/tH, where d7V./dt ~ iV,/[log(2/i?ic)tr] 
is the scattering rate into loss-cone orbits that take an SBH 
into the MBH, ~ w^/ (log AG^pm.) is the relaxation time 
(assuming all the mass is in SBHs), i?ic is the loss-cone open- 
ing angle and A ^ 0.4A'^, is the Coulomb cutoff (e.g., Syer 
& Ulmer [T999h . This constraint yields the relation (solved 
numerically). 



maxA'',(< Ti,) 



21og(2^r,/r,) ( mV P{r,) 



31og(0.4maxA^,) \m, J tn 



(1) 



where rg is the Schwarzschild radius of the MBH and P{ri,) 
is the orbital period at r*. 

3. THE DIRECT EXCHANGE CROSS SECTION 

We estimate the efficiency of the capture mechanism with 
the assumption that the angular momentum distribution of the 
young stars has been efficiently randomized. For simplicity, 
orbital periods and periapses are considered Keplerian (this is 
reasonable inside 2.5 pc, where the enclosed mass is < to) 
and the typical orbital energy of a star is represented by its 
virial energy (equivalent to assuming the star was initially on 
a circular orbit). If the star is scattered to an eccentric orbit, it 
will pass by the MBH on a hyperb olic orbit relative to it, with 
energy (Alexander & Livio l2001l) 
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where ao is the orbital semi-major axis (SMA) and /io"^ is the 
total mass enclosed within ao ( the tilde symbol denotes quan- 
tities in dimensionless units G=Mi,=Ri,=l; in these units, the 
stellar binding energy is ^ 1). 

A three-body encounter (Heggie'197;S) is characterized by 
a dimensionless relative velocity = —Eq/E, between the 
incoming star at infinity and the binary barycenter, where 
E, — ~ffifri,/2a, is the binding energy of an MBH-SBH 
"binary" with SMA a,. When Voo > 1 (a fast encounter with a 
soft binary), the incoming star carries enough energy to ionize 
(disrupt) the binary. When Voo < 1 (a slow encounter with a 
hard binary), the only possible outcomes are exchange, where 
the incoming star ejects one of the binary members and re- 
places it, or scattering, where the incoming star remains un- 
bound. Typical three-body encounters between an ^ 10 Mq 
star and a SBH in the GC occur below the ionization threshold 
(Table[0. Note that the actual 2-body interaction between the 
star and the SBH is analogous to the high velocity exchanges 
studied by Hut (1983). 

An exchange can proceed via two channels: direct or reso- 
nant. It is direct when the star passes within the SBH capture 
radius Fc ^ [(1 + TO,)/m]a, and ejects it. In this case, the 

— 1/3 

three-body analogy is justified because <IC , the mean 
distance between stars. It is resonant when a transient three- 
body bound system forms and persists for many orbits over 
a volume of radius ^ a,, until one of the masses is ejected 
(Hut 1993). In this case, the isolated three-body analogy is 

— 1/3 

no longer valid because a, ^71* . The direct exchange 
cross-section, which is the relevant channel here, is smaller 
by a factor of [(1 + to,)/to]'^ than the resonant one (Heggie, 
Hut & McMillan 1996) and was therefore neglected in past 
works on isolated three-body encounters. We derive it using 
the approximate equality between the binding energy of the 
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original and exchanged binary (Heggie et al. 119961), and we 
calibrate it by detailed numerical 3-body simulations (Fig. 
In the limit m ^ m,,l and for 1/2 < m, < 2, the direct 
exchange cross-section averaged over all orbital angles for a 
thermal distribution of binaries with SMA a, is 
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where the first factor is the gravitational focusing term, the 
second is the geometrical term irr^, and the third is the phase 
space factors required by detailed balance arguments (Heg- 
gie et al. 1996). The best fit parameters obtained from the 
numerical simulations are and i?~0.7. 

Equation (|3} describes an isolated system in which the bi- 
nary has a specified SMA and in which the captured star is 
initially unbound at infinity. Capture by three-body exchange 
in the CG differs in that the stars are scattered from a fi- 
nite distance and are captured by a distribution of SBHs en- 
closed within 7\. Equation (|3} is then applied by replacing 



a, (a.) 

replacing 



1.5r,t, the mean SMA in a r ^/"^ cusp, and by 
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whereuQ — —{iiQni/ao)/{2E,). The ionization threshold re- 
mains at = 1, but now ionization implies the ejection of 
the SBH out to a distance ~ oq, where its orbit is randomized, 
and not to infinity. 

The number of captured stars orbiting the MBH is the prod- 
uct of the capture probability per passage, A^.E, the typical 
lifetime of a star after capture, t*/2, and the incoming flux 
of young stars, [rii^qi,/ Pi^)d?r / (TrrJ), where Pq is the or- 
bital period at ao and = 1 — (1 — r*/ao)^ is the fraction of 
stars on orbits that cross inside r,t. The steady state number 
of captured young stars within r^, is then 



■ dao . 



(5) 
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where I'i, is the stellar mass density per unit SMA.^ The inte- 
gration runs from the smallest SMA where young stars can be 
formed (here assumed to be ai = 0.44 pc, where a star with 
the virial energy is unbound to the MBH) up to a2 — 2.5 pc 
where the Keplerian approximation is marginally valid. Stars 
can also be captured from orbits with aq > 02 up to a maximal 
SMA, where Voo ~ 1- However, the contribution of such stars 
to the total is small because of the longer Pq and smaller q^,. 
Note that since Voo increases with a, and ao, there exists a 
maximal radius where stars can be captured, < 1 pc. Be- 
yond this limit three-body interactions result in ionization, not 
capture. 

4. EXCHANGE CAPTURE IN THE GALACTIC CENTER 

Table Q shows the number of captured stars, assuming that 
SBHs born in the central 5 pc are now concentrated inside 
r, = 0.7 pc. This choice of r, (e.g. Morris 1993; Miralda- 
Escude & Gould 2000) implies that the central concentra- 
tion of SBHs is near the drain limit (Eq. ^ and that SBHs 



^ The spatial distribution p*(r) : 
sponds to Ui, (a) = 3it2^°' l3{a + l, 



: Cr in a Keplerian potential corre- 
'3/2)Ca^-°' fSch6deletal. l200l . 




Fig. 1. — Direct exchange cross-section below the ionization threshold 
(multiplied by v'^ /iral ) for m, = 2/3, 1 and fh = 16, 64 as a function 
of V. Points with error bars are Monte Carlo estimates of the cross section 
(using the SIGMA3 program [McMillan & Hut 1996] of the Starlab soft- 
ware package by Hut, McMillan, Makino & Portegies Zwart; available at 
http://www.ids.ias.edur starlab). Thick hues are the best-fit analytic approx- 
imation (Eq.0. hor D > 1, binary ionization competes with direct exchange, 
and the cross section (shown here only for m : m, : Af* = 16 : 1 : 1) rapidly 
falls off. 



Table 1 

Three-body exchange in the Galactic center 



Parameter 



Extent of young stars, r* (pc) 
0.04 0.10 



Number of observed young stars, A'^* 


~10 


~40 


Number of captured young stars in r*, A'^* " 


2.4 


7.7 


Number of SBHs in r^,, N.'' 


4.2x10^ 


1.3x10"* 


Drain limit, max A*'. 


4.7x10^ 


1.7x10"* 


Enclosed SBH / stellar mass'^ ratio in r* 


2.1 


1.5 


Mean'* velocity at infinity, (-fioo) 


0.2 


0.4 


Median'' initial SMA, ao (pc) 


0.9 


0.9 



" For 2.3 X 10 young stars (Mq = 3-15 Mq) inside 2.5 pc. 
For 1.5 X 10^ SBHs of 7 Mq within r, = 0.7 pc. 
Relative to the estimated enclosed stellar mass (Genzel et al. l200: 
dA'^* -weighted mean between ai and a2. 

^ Half of the captured stars originate between median ao and 02 



dominate the central mass density, as expected in a dynam- 
ically evolved system (Freitag 2003). The predicted mean 
number of captured stars is ^ 25% of those observed. The 
captured stars are expected to lie mostly at the low-mass and 
low-luminosity end of the observed range because lower-mass 
stars are more numerous and live longer 

In addition to the fact that capture by exchange can account 
for at least some of the young stars, this mechanism can nat- 
urally explain observed trends in their orbital properties. Be- 
cause Tc falls with the distance to the MBH, there exists a 
minimal distance rmin ^ 0.8mm,' /(I + m,) where the cap- 
ture radius equals the tidal radius of the SBH (Alexander & 
Kumar 2001). Beyond this limit, an exchange is no longer 
possible since the star is disrupted. Near disruption at Tmin, 
the tidal interaction absorbs orbital energy of the order of the 
stellar binding energy, so that the captured star loses most of 
its kinetic energy in the encounter, and r„iin then becomes the 
apoapse of the new orbit. In the GC, r,„in 0.01 pc, which is 
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consistent with the observed lower bound on the apoapse of 
the young stars (the local relaxation time is longer than i*, so 
they are expected to remain on their original orbits). We note 
that the extraction of orbital energy can increase the capture 
efficiency. A rough estimate indicates a factor of ^ 2 increase 
in the cross-section for captures occurring near r^nin, where 
the tidal interaction is strong. This effect was not included in 
the results quoted in Table Q. 

The numerical simulations also show that the exchanged bi- 
naries have higher than random eccentricities, even when the 
SBHs have a random (thermal) distribution of eccentricities. 
The bias becomes substantial the faster (voo 1) and more 
massive (m, < 1) the incoming star is. This property of the 
exchange mechanism agrees with the observed trend. 

Lower mass main-sequence stars in the range — 1- 
3 Mq (-G2V-A0V) can be similarly captured by NSs. Dy- 
namical simulations of the GC indicate that the central num- 
ber density of NSs is 1/3 that of the SBHs (M. Freitag 
2003, private communication). We find that there are ~ 35 
lower-mass stars captured in the inner 0.04 pc and ^110 cap- 
tured in the inner 0.1 pc. 

5. DISCUSSION AND SUMMARY 

The efficiency of a "billiard ball" recoil depends strongly 
on the mass ratio of the colliding objects. The dynamical 
evolution of a stellar system around an MBH naturally pro- 
vides a dense concentration of targets whose masses are well 
matched for stopping and capturing unbound young stars like 
those observed very near the MBH in the GC. The attempt to 
predict the number of captured stars is limited by the uncer- 
tainty in the mass and number distribution of the SBHs and 
in the number and orbital properties of young stars (spectral 
types '^BQV-OSV) in the inner few parsecs of the GC. We 
show that under favorable conditions, capture can account for 

25% of the observed young stars. 

Additional effects that were not taken into account here may 
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increase the capture efficiency. We considered only capture 
by a single direct exchange in the point-mass approximation. 
However, a few weaker interactions may also lead gradually 
to capture during the star's lifetime. Further study is needed 
to estimate the contribution of multiple scatterings to the cap- 
ture cross section. Tidal energy extraction in captures near 
the tidal limit also increases the capture cross section there. A 
possible outcome of internal mixing by a strong tidal interac- 
tion is an extended main-sequence lifetime and higher lumi- 
nosity (Maeder & Meynet 2000). This would further increase 
the number of captured luminous stars (Eq. |5j- 

In the context of direct exchange capture, the stars with the 
smallest apoapse offer an opportunity to study the long term 
effects of a strong tidal interaction. Interestingly, S2, the star 
with the smallest apoapse, is also the brightest (Ghez et al. 
2003; Schodel et al. 2003). The orbital eccentricity distri- 
bution reflects the mass ratio between the star and the SBH. 
Thus, a spectral determination of the masses of the captured 
stars together with a statistical analysis of their eccentricities 
can probe the poorly known mass function of SBHs. The to- 
tal number of captured stars over the lifetime of the Galaxy 
(assuming steady state), 2Ni,tH /ti, ~ 3 x 10^ (for 40 stars 
of 10 Mq in 0.1 pc), is of the same order as the number of 
SBHs in the inner parsec. Thus, if the dynamical friction 
timescale at (ao) (where the SBHs are ejected to) is not much 
smaller than the age of the Galaxy, the continual replacement 
of SBHs by NS progenitors and o f NSs by white dwarf pro- 
genitors (Fryer & Kalogera 1200 II) may regulate the build-up 
of the dense cusp of compact steflar remnants. 



We thank P. Hut for his generous help and R. Genzel, B. 
Hansen, D. Figer and M. Freitag for useful discussions. T.A. 
is supported by ISF grant 295/02-1, Minerva grant 8484 and 
a New Faculty grant by Sir H. Djangoly, CBE, of London. 



Kim, S. S. & Morris, M. 2003, ApJ, 597, 312 

Krabbe, A. et al. 1995, ApJ, 447, L95 

Maeder, A. & Meynet, G. 2000, ARA&A, 38, 143 

McCUntock, J. E. & Remillard, R. A. 2004, in Compact Stellar 
X-ray Sources, in press, ed. W. H. G. Lewin & M. van der 
Klis (Cambridge: Cambridge University Press), available at 
http://xxx.lanl.gov/abs/astro-ph/0306213 

McMillan, S. L. W. & Hut, R 1996, ApJ, 467, 348 

Mezger, P G., Zylka, R., Philipp, S., & Launhardt, R. 1999, A&A, 348, 457 
Miralda-Escude, J. & Gould, A. 2000, ApJ, 545, 847 
Morris, M. 1993, ApJ, 408, 496 

Mouawad, N., Eckart, A., Pfalzner, S., Schodel, R., Moultaka, 
J., & Spurzem, R. 2004, A&A submitted; available at 
http://xxx.lanl.gov/abs/astro-ph/0402338 

Portegies Zwart, S. R, Makino, J., McMillan, S. L. W., & Hut, P 2002, ApJ, 
565, 265 

Portegies Zwart, S. R, McMillan, S. L. W., & Gerhard, O. 2003, ApJ, 593, 
352 

Schodel, R., Ott, T., Genzel, R., Eckart, A., Mouawad, N., & Alexander, T. 

2003, ApJ, 596, I0I5 
Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 
Syer, D. & Ulmer, A. 1999, MNRAS, 306, 35 
Vollmer, B. & Duschl, W. J. 2001, A&A, 377, 1016 
Zhao, H., Haehnelt, M. G., & Rees, M. J. 2002, New Astronomy, 7, 385 



